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Abstract 

The application of probability density function (pdf) turbulence models is addressed in this work. For 
the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional CFD 
flow solver with the Monte Carlo simulation of the pdf evolution equation has been developed. The algorithm 
has been validated using experimental data for a heated turbulent plane jet. The study of H 2 -F 2 diffusion 
flames has been carried out using this algorithm. Numerical results compared favorably with experimental 
data. The computations show that the flame center shifts as the equivalence ratio changes, and that for the 
same equivalence ratio, similarity solutions for flames exist. 


Introduction 


It is the consensus of the combustion profession 
that the prediction of chemical reaction rate (the 
source term in a species conservation equation) is 
poor if a conventional turbulence model is used. The 
main difficulty lies in the fact that the reaction rate 
is highly non-linear, and the use of averaged tem- 
perature, pressure and density produces excessively 
large errors. Moment closure models for the source 
terms have attained only, limited success because the 
assumptions for such models to be valid often can 
not be satisfied. The probability density function 
(pdf) method seems to be the only alternative at the 
present time that uses local instantaneous values of 
the temperature, density, etc., in predicting the chem- 
ical reaction rate, and thus is the only viable approach 
for more accurate turbulent combustion calculations. 
Two main lines are being followed in pdf methods: 
one uses an assumed shape for the pdf, the second 
solves a pdf evolution equation; the present paper ad- 
dresses only the latter. There has been significant 
progress in the study of pdf turbulence models in low 
speed flows in the past decades. These developments 
were summarized in Refs. [1], [2] and [3]. In spite 
of this progress, pdf turbulence modeling remains a 
nascent discipline with many unresolved issues. 

The fact that the pdf equation has a very large 
dimensionality renders the use of finite difference 
schemes extremely demanding on computer memories 
and thus impractical, if not entirely impossible. A log- 
ical alternative is the Monte Carlo scheme, which has 
been used extensively in statistical physics. The evo- 
lution equations for the joint pdf of the velocity and 
species mass fraction have been successfully solved 
using Monte Carlo schemes, see, e.g., Pope[l]. How- 
ever, since CFD has reached a certain degree of ma- 
turity as well as acceptance, it seems, at least from 
the standpoint of practical applications, that the use 
of a combined CFD and Monte Carlo scheme is more 
beneficial. Therefore, in the present study a scheme 
is chosen that uses a conventional CFD algorithm to 
solve the Navier-Stokes equations and provide flow- 
field properties such as velocity, pressure, etc., and the 
chemical reactions are calculated by using a Monte 


Carlo scheme to solve a pdf evolution equation. 

The combined CFD-pdf solver has been devel- 
oped recently and validated using non-reacting flow 
data for a heated turbulent plane jet. The algorithm 
has been further tested in the numerical study of an 
Hi-F? diffusion flame. This diffusion flame was stud- 
ied experimentally by Mungal k Dimotakis[4] and 
Hermanson k Dimotakisfo]. The numerical results 
from the present study are compared with these ex- 
perimental data. 


Theory 

Governing equations for reacting flows. 

Flows with chemical reaction are governed by 
the continuity equation, momentum equations, en- 
ergy equation, and species transport equations: 


d t p + dipui = 0 

pd t u, -t pujdjUi ~ -dip + pdjTij 

pd%h 4- pujdj - d t p — Ujdjp 

— —djqj 4- Q 

pdtYk 4- pujdj Yk = pdj(DdjYk) 4- w k 
k = 1, 2, • - - , A\ 


where u, is the velocity, V* is the mass fraction, and 
Wk is the chemical source term. (In addition to these 
equations, one also needs the equation of state.) For 
turbulent flows, we substitute 


Ui - u, + u ' , 
Yi = Yi 4- Yf, 


(2) 


into the above equation and take an ensemble aver- 
age. In the process of averaging, new unknown quan- 
tities in the form of correlations appear, e.g., uju'-, 

uJ Y*, etc. These quantities can all be modeled using 
conventional turbulence models, such as two equation 
models or second order closure models. A problem 
unique to flows with chemical reaction is the average 
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of the chemical source term, puT*, The major diffi- 
culty lies in the fact that Wk is a highly nonlinear 
(usually exponential) function of the temperature. It 
is well know that the use of averaged temperature, T, 
in evaluating pwi can cause egregious errors. One may 
consider the effect of the temperature fluctuation, T\ 
by applying a moment closure model to the term ~pwl\ 
however, such a closure model results in an infinite se- 
ries that converges only when T a T and T* T. In 
many combustion problems, these two conditions are 
violated: in fact, we often have T a ^ T and T* T 
instead. In view of the above, the prospect of an ac- 
curate prediction of turbulent combustion using con- 
ventional turbulence models seems dismal. The above 
fact motivated the use of pdf methods. 

Pdf evolution equation 

Given a set of m random vari- 
ables, and the joint probability den- 
sity function, P^i. V> 2 * * * • , y, z % t), the mean of 

any random function, can 

be calculated as 

f(x,y,z,t) = J ■ J f P difri dfa - ■ ■ (3) 

For simple 'Hows, one could assume the shape of 
the pdf and compute the source term, ptZT, based on 
the assumed pdf using the above integral. For more 
general problems, one needs to solve a pdf evolution 
equation for a more accurate pdf distribution. In the 
latter case, the evaluation of the source term, pu7jt, is 
no longer necessary because the mean values of the 
temperature, species mass fractions, which we previ- 
ously would obtain by solving the transport equations 
(Eq. 1), can now be evaluated directly from the pdf 
using the above integral. 

The pdf evolution equation can be derived from 
the transport equations (Eq. 1) in many different 
ways, using a Dirac delta function, a characteristic 
function, or a characteristic functional. [1-3] The evo- 
lution equation of a single point probability density 
function of scalar random variables tpi , * * • , 0 m can be 
written as 

m 

pd,P + pv a d a P + p ,'Pm)P] 

1 = 1 

= -d a (p < v'Mx) = r/, t > P) 

m m 

-pEE^« «<#(**(*) = t h>p) (4) 

.= 1>=1 

where the terms represent the rate of time change, 
mean convection, chemical reaction, turbulent con- 
vection, and molecular mixing, respectively; P is the 
density- weighted joint pdf: 

P = PP/P . 

(ij is the scalar dissipation: 

Uj = Dd Q <t>id Q 4> h 


(where D is the diffusion coefficient), and < rjy > 
denotes the mathematical expectation of a random 
function x conditioned upon y. 

The left hand side of eq. (4) can be evaluated 
exactly and requires no modeling; the right hand side 
terms contain the conditional expectation of the ve- 
locity fluctuation and the conditional expectation of 
the scalar dissipation, which are new unknowns and 
require modeling. 

Closure models for pdf cnuation 

The first term that needs modeling is the tur- 
bulent convection term, for which one can use the 
following gradient model: 

-<t>”crhhfc > P — D t d a P t 

where D t is the turbulent diffusion coefficient, which 
is set to be equal to the eddy viscosity, i.e., the tur- 
bulent Schmidt number is equal to one. 

The next term needs modeling is the molecular 
mixing term. A coalescence/dispersion model is used 
for this term, which has the following general form[6]: 

N /v 

-^EE4,(< ( .i \+k>P) 

;=i i = i 

j J *"Wd+" - p ' 

= M{P) 

where T is the transition probability. A new mix- 
ing model continuous in time is recently developed by 
Hsu and Chen[7]. For more detail description of the 
molecular mixing models, see Refs. 6 and 7. 

The modeled pdf equation is then 

m 

pd,P + pv a d a p + p . • • • , iPm)P} 

1=1 

=dj(D t d a P) + M (P), (5) 

which can be solved using a Monte Carlo simulation. 

Numerical Methods 

In the present study, a combined finite 
difference— Monte Carlo solver is developed. For the 
velocity field, the Navier-Stokes equations and a k — i 
turbulence model are solved using a finite difference 
method; for scalar variables such as the temperature, 
mass fractions, etc., the pdf evolution equation is 
solved using a Monte Carlo scheme. The CFD flow 
solver provides the Monte Carlo solver with the mean 
velocity and a turbulence time scale r, where r = k/c, 
and the Monte Carlo solver provides the mean flow 
solver with the density, p. 

The solution of the Navier-Stokes equations 
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Since the primary concern of the present work is 
pdf modeling, we choose the simplest possible solver 
for the N-S equations. The equation is transformed 
into a general coordinate system using the following 
general relation: 

x = { 

y ~ s/(£. m 


d v = —d„ 

yn 

The transformed momentum equation for steady 
flows is 

ud(ii + (v - ut/()—d q u = ~^d(p + —d n (v T d n u), 
Vij P Vr) 

where the viscous term involves gradient in the {- 
direction is neglected since we are only interested in 
shear flows in the present study. 

First order upwind difference is used in the 
{-direction so that a marching scheme for steady 
parabolic flows can be used. To ensure stability, a 
flux splitting scheme is used in the redirection. 


(3) chemical reaction: 

P n+l =(1 + A tR)P" 

In the present study, only steady flows are considered, 
so dt is replaced by d{ and a marching scheme in the 
{-direction is employed. 

In a Monte Carlo simulation, the continuous pdf 
is replaced by N x M delta functions, 

n— I 

x 6 (xp 2 - 4 n> (0) • ■ • Wm-rfSftO). 

where each product of the m delta functions repre- 
sents one event of an ensemble of N sample events. 
An event can be thought of as a fluid particle, and 
the evolution of P * entails the movement of the par- 
ticles in the physical space as well as the phase space 
(^— space). The movement of the particles is. of 
course, governed by the pdf evolution equation, and 
is simulated in the following three steps. 


The solution of the pdf equation 

A fractional step Monte Carlo method, as will 
be described in the following, is used to solve the pdf 
equation. We first discretize the time derivative in eq. 
(5) using finite difference: d t P = (P n+l - P n )jAt , 
then the pdf evolution equation (eq. 5) can be written 
as 

m 

p n+i - {1 - At v a d a - At d a ,wi(4> i, ■ ■ ■ . i p m ) 

i=i 

+Atd a D,d a + At M } P n 


Step I: convection 

Replace P by P* in the pdf evolution equation 
and transform it into the { — r; coordinate system, the 
convection process for a steady flow can be written as 

pud{ P* — — p{v — uy^—drjP' 4 d n (pD t d r7 P m ) i 

Vt) y n 

Discretizing the above equation using a finite differ- 
ence scheme, we can write 


Using approximate factorization, the above equation 
is recast as 


p:, = °Pi- 


ij+i 


■3P’. U 


rPUi-i 


p n+1 =(1 - At yx 

i—l 

x ( 1 + A*,\/)(l - A tpvd a -r A td Q D t d Q ) P n 

+0(At) 

= (1 + AfC)(l + AtM)(l 4- A tR)P n 4- O(At), 

where C denotes the convection operator, m the 
molecular mixing, and Ft chemical reactions. With 
the above expression, we can three processes consec- 
utively: 

(1) convection: 

P m = (I + A tC)P n , 

(2) molecular mixing: 

P** = (1 + A t\f)P\ 


(Here again we used a one sided difference so that a 
marching scheme could be used.) The above equation 
states that if we divide the flow field into cells, then 
the pdf at point (2,7) can be written as a linear com- 
bination of the pdf at neighboring points, To simulate 
this process with a Monte Carlo scheme, we move the 
sample particles between cells according to the above 
equation. For instance, the particles of cell (1,7) will 
be obtained by choosing randomly aN particles form 
cell (*- 1,7 4- 1), PN particles from cell (i- 1,7), and 
jN particles from cell (1 — 1,7 — 1). In order for the 
total number of particles not to change, we require 
that a 4- 0 4- 7 — 1. 

It is worth noting that this method can easily 
be extended to elliptic flows and applied to general 
curvilinear coordinate systems. 

Step 2: molecular mixing. 
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The molecular mixing process is simulated bv the 
following binary interaction model: 

pud^P 

=^r / / - p 

where T is the transition probability. By assigning 
various functions to T, we would have different mix- 
ing models. In the present study, the modified Curl 
model by Janicka et al. [8] is used. In this model, the 
transition probability is given as 


will acquire new composition at the next moment, and 
the changing rate is 


i- 1,2, - , Af, 

where w t is the reaction rate for that specific sample 
point. If we regard u; t as the convection velocity of 
the particle in the phase space, the it is clear that 
the above ODE represent the movement of a sample 
particle in the phase space. 

The various chemistry models considered in this 
study are presented in the following section. 




for tl>'<il>< 
or ip” < ip < 
0 otherwise. 


with Cd = 6.0. 

In the Monte Carlo simulation, the above model 
is realized in the following manner: Divide the flow 
domain into small cells, each containing N sample 
particles. Given a small time interval At and a tur- 
bulent time scale r, select randomly /V mr pairs of par- 
ticles, where 

N mx = Q.b^-N, 

Ct 

and let a pair, say, m and n. mix as follows 

6 n (t + At) = A4> m {t) + (1 - A)4> n {i) 

+ At) = A<fi n (t) + (1 - .4)0 m (t) 


Chemistry Models for /A - F-, Reaction 
Finite rate chemistry . 

The reaction of H 2 + F 2 = 211 F can be repre- 
sented by the following two step chain reactions^]: 

H 2 + F — - kl HF + H 

H + F 2 — .*» HF+F 


with 

= 2-6 x 10 12 T° 5 exp , 


i 2 = 3x 10 9 T ls exp 



where T is in K , R in cal mol 1 K~ l t and k in cm 3 
mo/ -1 sec -1 . 

The above expressions only provide the forward 
reaction rate; the backward reaction rate can be cal- 
culated using the equilibrium coefficient: 


where A = 0.5£, with £ a random variable uniformly 
distributed on the interval [0.1). The remaining N - 
2A r mx particles remain unchanged: 

+ AO = o n (0 

The turbulent time scale is supplied by the finite 
difference solution of a £ - c turbulence model: r = 
k/c. 

This model admits the non-physical jump con- 
dition and does not produce the correct long time 
behavior for decay problems in homogeneous turbu- 
lence. A continuous model that predicts the correct 
long time behavior for turbulence decay problem has 
been introduced by Hsu and Chen [7]. But as shown 
in Ref. [7], for practical combustion problems where 
the long time statistical behavior is not crucial, the 
modified Curl model gives acceptable results. In the 
present study, both the modified Curl model and the 
continuous model have been used. 

Step 3: chemical reaction 

Chemical reaction is represented in the Monte 
Carlo simulation by the movement of sample particles 
in the phase space due to reaction. A sample point 
with a given composition {oi,4 3t - • • ,4 m } at time t 


f i 

and the equilibrium coefficient can be evaluated using 
Gibbs free energy. Polynomials for Gibbs free energy 
for various species can be found in Ref. [8). A nu- 
merical experiment shows that the backward reaction 
rates are much smaller than that of the forward re- 
action for the two reaction steps described above and 
can be neglected. 

Let C denote the mole concentration; for each 
sample particle, we need to solve the following set of 
ODEs: 

dC H, _ , „ r 


dt 

dCr 7 
dt 
dC H 
dt 
d Cf 


— — k? Off Cf 7 — ^2 


— — <-^2 — 


dt 
dCffF 
dt 


— U)\ — U>2 = CJ4 


— — U>\ — u>2 — - tJs 


4 A V- ,.0 

1=1 


where h is the .enthalpy and h is the heat of forma- 
tion. 



Using n to denote the time level, we found that 
the following semi-implicit scheme seems to be most 
stable: 


a few hundred degrees Kelvin. Under these circum- 
stances, the variation of specific heat C T is negligible, 
and the temperature can be calculated as 


£"1 + 1 _ Cff, 

1 + A <jb(C£ 

C n 

C n+1 = UFl 

F> 1 + At kiC n H 

, = Q + A tkrC^Cl 
H 1 + A t k-iC n H 

_ Q + A^,Q,C^ 

F l + Alt.C^ 

C"hf =C n HF + (kiC^Q. + k 2 Cr H C" Fi ) 

The energy equation is solved explicitly. To start the 
chain reaction 0.00002 mole fraction of F was released 
into the flowfield initially. 

Using the above equations to calculate the reac- 
tion in a laminar premixed flame, we found that the 
time required for this reaction is the order of one mi- 
crosecond (Figs. 1 and 2). The turbulence time scale 
in the calculation, r = Jc/e, is the order of one sec- 
ond. Considering this large difference of time scales, 
it is impractical to use finite rate chemistry in this 
calculation. 

Fast reaction 

Since the chemical reaction is very fast compared 
to the flowfield development and the reverse reactions 
are negligible, we chose to use the following complete 
irreversible reaction mechanism. We assume that af- 
ter molecular mixing, complete reaction is achieved 
within each particle during the time interval of one 
marching step in the calculation. The reaction is cal- 
culated using the following equations. 

For Ch 7 > Cp 3 : 

U H 7 - 

eg 1 = 0 

S~*n i 

— ^HF + 

For 

+ 1 , r, r-m 

The energy equation can be written as 

A* = -*?„, A C„ F {^f). 

where Af// F if the molecular weight of HF. 

Since the concentrations of reactants are low, the 
heat release is low, and the temperature rise is within 


rpn + l _ 'j-n _j_ 


Ah 


The error involved in this approximation is less than 

5%. 


Results and Discussions 
Code validation. 

Before applied to the hydrogen-fluorine diffusion 
flames, the computer algorithm is first vaJidated us- 
ing experimental data for a heated turbulent jet. The 
non-reacting flow data serves as a check for the con- 
vection and molecular mixing process in the pdf solver 
as well as the k — c model in the N.-S. solver. 

Extensive experimental results for turbulent 
plane jet have been reported by many authors. Mea- 
surements for mean velocity field in a turbulent jet 
were first reported in the 1930’s[9]; turbulent shear 
stress measurements had been reported more re- 
cently [10, 1 Ij . To determine the effect of turbulence 
on mixing, the temperature field of a heated turbu- 
lent jet had been studied by several authors. The tur- 
bulent jet has a slightly higher temperature than the 
ambient. Measurements of both the the mean tem- 
perature and the rms of the temperature fluctuations 
were given[12-17]. 

In the present study of this non-reacting flow 
case, the temperature field is treated as a conserved 
scalar and is simulated by the pdf of the tempera- 
ture; the velocity field and turbulent shear stress are 
obtained by solving the N.-S. equation and a k — i 
turbulence model. 

In the finite difference solution of the flowfield. 
41 grid points are used across half of the jet width. A 
symmetry boundary condition is used at the jet cen- 
terline. Fig. 3 shows the comparison of the present 
solution of the mean velocity field with the experimen- 
tal data, and Fig. 4 presents the numerical solution 
of the turbulent shear stress as compared to the ex- 
perimental data. Good agreements between numer- 
ical results and experimental data are observed for 
both the mean velocity field and the turbulent shear 
stress. The turbulent shear stress is calculated from 
— < tiV >= u T du/dy , where v r = Cp/^/e. A good 
prediction of the turbulent shear stress ensures that 
the turbulent time scale, r = k/e, supplied to the 
Monte Carlo simulation is correct. 

For the Monte Carlo simulation of the temper- 
ature field, two sample sizes of 1000 and 1500 sam- 
ple particles per cell are used. The predicted mean 
temperature and the root mean square (rms) of the 
temperature variation are presented in Figs. 5 and 
6. The results show that both calculations produce 
fairly good comparisons with the experimental data, 
which means a sample size of 1000 particles per cell 
is large enough.* 
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The above validation lends credibility to the N.S. 
solver as well as the Monte Carlo solver developed in 
the present study 

Hi - A diffusion flames 

The flow conditions for the Hi — A diffusion 
flames are set according to an experiment performed 
by Hermanson and Dimotakis (1989). The flame con- 
sists two streams. The upper stream contains A and 
Fi, the flow velocity is U\ = 22 m/s; the lower stream 
contains of Ni and of Hi, with velocity U 1 = 8.8 m/s. 
In the present study, 6 cases involving various per- 
centage of H? and A in the upper and lower stream 
are considered; the conditions are listed in Table l t 
where A T*# is the difference between the adiabatic 
flame temperature and the free stream temperature. 


ox 

No. 

equivalence 

ratio 

lower stream, 
mole fraction 
of H 7 

upper stream, 
mole fraction 
of A 

AT.* 

i 

l 

0.02 

0.02 

186 

2 

1 

0.04 

0.04 

368 

3 

1 

0.06 

0.06 

554 

4 

1/4 

0.01 

0.04 

151 

5 

1/4 

0.02 

0.08 

302 

6 

1/4 

0.04 

0.16 

ESI 


Table I. Initial conditions of the flames calcu- 
lated in the present study. 


Fig. 7 shows the calculated temperature rises due 
to combustion for cases 2 and 5 of Table 1 and the cor- 
responding experimental data. In the figure, 67 is the 
shear layer thickness determined by 1% of the tem- 
perature rise, AT is the actual temperature rise due 
to combustion. ( the two streams have the same tem- 
peratures initially,) and A T a dj is the adiabatic flame 
temperature assuming complete reaction. The solu- 
tion for case 2 agrees fairly well with the experimental 
data for the same case, while the mean temperature 
is slightly over predicted for case 5. The computation 
shows the flame center shifts as a result of change in 
equivalence ratio, which is consistent with the exper- 
imental results. 

The rms of the temperature variance for cases 2 
and 5 are given in Fig. 8. The results show that in 
a diffusion flame, the temperature variance has two 
peaks, and the highest values do not coincide with 
the maximum values of the temperature distribution. 
The shifting of the flame due to the change of equiv- 
alence ratio can also be observed from this figure. 

The mass fractions of Hi, Fi and HF for cases 
2 and 5 are presented in Figs. 9 and 10, respectively. 
One can see that although a complete reaction chem- 
istry model was used, with the pdf method, results 
simular to that of a finite rate computation are pro- 
duced, which is one of the many advantages of the 
pdf method. 

One experimentally established fact is that for 
the same equivalence ratio, with various mole con- 
centrations of fuel and oxidizer in the flowfield, the 


normalized temperature stays the same [4,5]. The 
present calculation confirmed this. Figs. 11 and 12 
are the mean temperatures and rms’s of the temper- 
ature variance for an equivalence ratio of one; three 
mole concentrations of fuel were considered. One can 
see that the curves coincide with each other. The 
same agreements are found for equivalence ratio 1/4 
(Figs. 13 and 14). 

The pdf distributions for flame temperature at 
the center and at the outer edge of the flame (for case 
2 of Table 1) are plotted in Fig. 15, where the x- 
axis denotes (T— < T >)/ cr , with < T >= T being 
the mean temperature and cr the rms. The pdf dis- 
tributions show that at the center of the flame, the 
temperature with the highest probability is not far 
from the adiabatic flame temperature, while at the 
outer edge of the flame, most of the time one would 
find a temperature close to that of the free stream 
temperature. Nonetheless, as a result of external and 
internal intermit tency, low temperature fluid does ex- 
ist at the center and high temperature fluid the outer 
edge. Pdf distributions for various species concentra- 
tion can also be obtained from the solution, but will 
not be presented here. 

Concluding Remarks 

A Monte Carlo solution algorithm for the pdf 
evolution equation has been developed and success- 
fully combined with a finite difference flow solver in 
the study of turbulent combustion. The algorithm 
was validated using turbulent mixing data from non- 
reacting flows. Turbulent diffusion flames of Hi- 
Fi were computed using the pdf method, and good 
agreements between numerical solution and experi- 
mental data were observed. The computation iden- 
tified the change of equivalence ratio as the cause of 
flame shift, and demonstrated that similarity solu- 
tion exists for flames with the same equivalence ra- 
tio. The present work showed that a grid depen- 
dent Monte Carlo scheme can be readily applied to 
a general curvilinear coordinate system: therefore, it 
is suitable for realistic reacting flow computations. 
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Figure 1. Temperature rise in a premixed lami- 
nar H 2 -F 2 flame. 



Figure 2. Changes of mass fractions in a pre- 
mixed laminar H ?-F? flame. H^, Fj, — - 

HF. 



y/y.su 0 


Figure 3. Mean velocity distribution in a 2D tur- 
bulent jet. o Ref (9], — Present calculation. 



Figure 4. 'IWbulent shear stress in a 2D turbu- 
lent jet. o Ref [10], A Ref [11], — Present calculation. 



Figure 5. Mean temperature in heated plane jet. 
A Ref [14], n Ref [17], o Ref [16] o Ref [15]. — pdf 
solution, 1500 particles per cell; - - - pdf solution. 
1000 particles per cell. 



Figure 6. RMS of temperature variance in heated 
plane jet. o Ref. [12], o Ref [17], A R*f (16], - pdf 

solution. 1500 particles per cell; pdf solution. 

1000 particles per cell. 
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Figure 7. Mean temperature rises in H 2 -F 2 
diffusion flames. Case 2 (refer to Table 1): o Ref [5], 
— present; Case 5: A Ref [5], present. 



y/<*r 

Figure 8. RMS of Temperature variance in // 2 - 
A diffusion flames. — equivalence ratio 1 (case 2, 
Table 1). equivalence ratio 1/4 (case 5). 



Figure 9. Mass fraction distributions in an 

F 2 flame, case 2, equivalence ratio 1. — // 2 , A 

- - HF. 



Figure 10. Mass fraction distributions in an tf 2 - 
Fi flame, case 5, equivalence ratio 1/4. — Hi, 

n rr 


Fi, — HF. 



Figure 11. Mean temperature rises in H?-F~ 
flames, equivalence ratio 1. Referring to Table 1: — 
case 1, o case 2, A case 3. 



Figure 12. RMS of temperature variance in H 7 - 
F 2 flames, equivalence ratio 1. Referring to Table 1; 
— case 1, o case 2, A case 3. 
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Figure 13. Mean temperature rises in H 2 -F 2 
flames, equivalence ratio 1/4. Referring to Table 1: 
— case 4, o case 5, A case 6. 



Figure 14. RMS of temperature variance in Hi~ 
F 2 flames, equivalence ratio 1/4. Referring to Table 
1: — case 4, o case 5, A case 6. 
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(T — <T > )/ 0 

Figure 15. Probability density function distribu- 
tions in an /fj-Fj flame (case 2). — pdf at the center 
of the flame, pdf at the out edge of the flame. 


86 







